Stable Equilibrium Based on Levy Statistics: Stochastic Collision Models Approach 
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We investigate equilibrium properties of two very different stochastic collision models: (i) the 
Rayleigh particle and (ii) the driven Maxwell gas. For both models the equilibrium velocity distri- 
bution is a Levy distribution, the Maxwell distribution being a special case. We show how these 
models are related to fractional kinetic equations. Our work demonstrates that a stable power-law 
equilibrium, which is independent of details of the underlying models, is a natural generalization of 
Maxwell's velocity distribution. 

PACS numbers: 

There is a strong analogy between the Gaussian cen- 
tral limit theorem (GCLT) and the relaxation to ther- 
mal equilibrium of the Boltzmann equation (y| and Ref. 
therein). However the GCLT is non-unique, which may 
imply that standard thermal equilibrium is non-unique. 
The Levy central limit theorem (LCLT) considers the 
case of summation of independent identically distributed 
random variables with an infinite variance Q. Hence 
following Montroll and Shlesinger it is natural to ask 
if generalized equilibrium concepts based on LCLT are 
meaningful Here we start answering this question 
using two very different types of collision models which 
still reveal the same type of equilibrium. We note that 
Levy statistics has many physical applications 0, how- 
ever its possible relation to generalized forms of equi- 
librium statistical mechanics is an open field of research. 
Recently, Bobylev and Cercignani |l( , investigated a non- 
linear Boltzmann equation with an infinite velocity vari- 
ance showing that the solution exists, and obtaining cer- 
tain bounds on it. In Q] the possibility of a relation be- 
tween solutions of the Boltzmann equation 5] and LCLT 
was briefly pointed out. 

Our goal in this Letter is to show that a new equilib- 
rium concept naturally emerges from old stochastic colli- 
sion models. Our models demonstrate that: (i) Levy ve- 
locity distributions serve as the natural generalization of 
the Maxwell velocity distribution, (ii) generalized power 
law equilibrium can be derived from kinetic models, there 
is no need to postulate a specific form of power law equi- 
librium, and (iii) the Levy equilibrium obtained here pos- 



sesses a certain domain of attraction, is unique, and does 
not depend on certain details of the underlying models. 

Model 1 We consider a one dimensional tracer particle 
with the mass M coupled with gas particles of mass m. 
The tracer particle velocity is Vm ■ At random times the 
tracer particle collides with gas particles whose velocity 
is denoted with v m . Collisions are elastic hence from 
conservation of momentum and energy = £i^m + 
&v m , where £i = fe, 6 = #r , e = m/M is the 



mass ratio and (V^) is the velocity of the tracer 
particle after (before) a collision event. The duration of 
the collision events is much shorter than any other time 
scale in the problem. The collisions occur at a uniform 
rate R. The probability density function (PDF) of the 
gas particle velocity is f(v m ). This PDF does not change 
during the collision process, indicating that re-collisions 
of the gas particles and the tracer particle are neglected. 

Many works considered this type of model, impos- 
ing the condition that the gas particles are distributed 
according to Maxwell's law, i.e f(v m ) is Maxwellian. 
Since we are now investigating possible generalizations of 
Maxwell's law we change this strategy and assume that 
f{vm) is non-Maxwellian. The goal is to see when and 
how the tracer particle reaches a universal equilibrium, 
which does not depend on the detailed shape of f(v m ). 

We now consider the equation of motion for the tracer 
particle velocity PDF W(VM,t) with initial conditions 
concentrated on Vm (0). Standard kinetic considerations 
yield the linear Boltzmann equation 



X 



W (V M ,t) = -RW (V M , T) + R dV 7 



M 



dv m W (V M ,t) f (v m ) x S (V M - £iV M - £, 2 v r , 



(1) 



where the delta function gives the constrain on energy 
and momentum conservation in collision events. Us- usual 
the first (second) term in Eq. (JJJ describes a tracer par- 



ticle leaving (entering) the velocity point Vm at time t. 
Eq. 0) contains convolution integrals in velocity space 
hence we consider now its Fourier transform (FT). Let 
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W(k, t) be the FT of the velocity PBF_W(V M ,t). Using 
Eq. QJ, the equation of motion for W (k,t) is a finite 
difference non-local equation 

W{k,t) = -RW(k,t) + RW(k^,t)f(k^), (2) 

where f (k) is the FT of f(v m ). The solution of the 
equation of motion Eq. JSJ) is obtained by iterations 



w(k,t) = j2 



{Rt) n exp (-Rt) 



n=0 



(3) 

This solution has a simple interpretation. The proba- 
bility that the tracer particle has collided n times with 
the gas particles is given according to the Poisson law 
P n (t) = £^ exp (-#*). Let W n {V M ) be the PDF of 
the tracer particle conditioned that the particle experi- 
ences n collision events. It can be shown that the FT of 
W n (V M ) is W n {k) = e ikv ^(°^U^J(k^-%). Thus 
Eq. @ is a sum over the probability of having n col- 
lision events in time interval (0, t) times the FT of the 
velocity PDF after exactly n collision event. 

In the long time limit W eq (k) = lim^oo W(k, t) an 
equilibrium is obtained from Eq. J2J). We notice that 
when Rt —> oo, P n (t) = (Rt) n exp(— Rt)Jn\ is peaked in 
the vicinity of (n) = Rt hence it is easy to see that 



W eq (k) = lim U^J(kQ-%) 



(4) 



In what follows we investigate properties of this equilib- 
rium. We note that similar equilibrium can be obtained 
also if the collision process is not described by the Pois- 
son law, any P n {t) which is peaked on n -* oo when 
t — > oo, with (nearly) zero support for finite values of n, 
will exhibit this behavior. 

We will consider the Rayleigh weak collision limit 
e — > 0. This limit is important since number of colli- 
sions needed for the tracer particle to reach an equilib- 
rium is very large. Hence in this case we may expect the 
emergence of a general equilibrium concept which is not 
sensitive to the precise details of the velocity PDF f(v m ) 
of the gas particles. In this limit we may also expect that 
in a statistical sense v m <C Vm, hence the assumption of 
a uniform collision rate is reasonable in this limit. 

We assume that statistical properties of gas particles 
velocities can be characterized with an energy scale T. 
Since T, m and v m are the only variables describing the 
gas particle we have 



/ («m) = 



1 



(5) 



We also assume that f(v m ) is an even function, as 
expected from symmetry. The dimensionless func- 
tion q(x) > satisfies a normalization condition 
<l(x)dx ~ 1, otherwise it is rather general. The scal- 
ing assumption made in Eq. (JSJ) is very natural, since the 
total energy of gas particles is nearly conserved. 



We first consider the case where moments of / (v m ) are 
finite. The second moment of the gas particle velocity is 
i^m) = m I^oo x2 Qi x )dx. Without loss of generality we 
set x 2 q(x)dx = 1. The scaling behavior Eq. (J3J) 
yields (v 2 ^) — (^) 52nj where the moments of q(x) are 
defined according to qm — x 2n q(x)dx. Thus the 
small k expansion of the gas particle characteristic func- 
tion is 



Tk 2 



?4 



Inserting Eq. 10 in Eq. (0} we obtain 



In [W eq (*)] 



T 
2m 



92 (e) k 2 



q A 



4! 



T 



g 4 (e) fc 4 +0(fc 6 
(7) 

where g n (e) = (2e)"/[(l + e) n - (1 - e) n ]. The interesting 
thing to notice is that in the limit e —> 0, the second term 
on the right hand side of Eq. J7J) is zero, thus q^ is an 
irrelevant parameter in the problem. In similar way one 
can show that all terms in the expansion containing q^ n 
with n > 1 vanish in the Rayleigh limit e — > 0. Thus 
using Eq. J7J 



lim In [W eq (k)] 



T 

2M' 



k 2 . 



(8) 



From Eq. (JHJ) it is easy to see that the Maxwell velocity 
PDF for the tracer particle M is obtained. Thus as ex- 
pected Maxwell's equilibrium is stable in the sense that 
for a large class of gas particle velocity PDFs Maxwell 
equilibrium is obtained. 

Now we assume that / (v m ) has a power law behavior, 
i.e., q(x) oc |x| _< - 1+Q ' when \x\ — » oo and < a < 2. For 
this case the gas particle characteristic function is 



/(*) 



1- 



r(l + a) 



T 



a/2 



- \k\ 



%3 



r(i + /3) 



T \ |8/2 

- Ifc^+odAsI") 
m J 

(9) 



where a < (3 < 2a. q a and qp are dimensionless numbers 
which depend of-course on q(x). Without loss of gener- 
ality we may set q a = 1. In Eq. @ we have used the 
assumption that f(v m ) is even. 

Inserting Eq. © in Eq. 10} we obtain a small k expan- 
sion In [W eq (fc)] . Taking the limit e — > one can show 
that the terms containing qp are much smaller than the 
leading term In [W eq (k)] oc — \k\ a . Thus qp, and in a 
similar way higher order coefficients, become the irrele- 
vant parameters of the problem. Thus we find that the 
tracer particle equilibrium characteristic function is 



W eq (k) ~ exp 



T 



aF(l-fa) \M 



a/2 



\k\ a 

A-a/2 



, (io) 
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FIG. 1: We show the FT of the equilibrium velocity distri- 
bution of the tracer particle. Numerically exact solutions of 
the problem are obtained using three long tailed gas particle 
velocity PDFs defined in text: case 1 squares, case 2 trian- 
gles, and case 3 diamonds. The tracer particle equilibrium is 
well approximated by the Levy distribution the solid curve; 



W eq (k) -exp (-2.211 
used T = 4.555, e : 



k |3/2 
TJs\ 
le - 5 and M 



For the numerical results we 
1. 



thus a Levy velocity distribution for the tracer particle 
is obtained. For a ^ 2 the equilibrium Eq. i|10|) depends 
on e, while for the Maxwell's case a = 2, the equilibrium 
is independent of the coupling constant e. Eq. (fTUl) im- 
plies that variance of the velocity diverges when a < 2. 
For the non-Maxwellian case the velocity distribution is 
characterized by the scale (T/M)"/ 2 which determines 
the width of the velocity distribution. 

The asymptotic behavior Eq. I|10|) is now demon- 
strated using numerical examples. We consider three 
types of gas particle velocity PDFs, for large values of 
\v m \ — ► co these PDFs exhibit / (v m ) oc |w m |~ 5 / 2 , namely 
a = 3/2. Case 1 f(v m ) = iV(l + 3 2/ V™/( 4T )|wm|) 5/2 , 
Case 2, f(v m ) = N 2 /[l + mvl l /{QA^T)fl i , where Nj. 
and N2 are normalization constants. Case 3, the gas par- 
ticle velocity PDF is a Levy PDF with index 3/2 whose 

< T n3/4 |fc|3/2 - 



T(5/2) 



FT is/(fc)=exp [-(X) 

According to our theory these power law velocity 
PDFs, yield a Levy equilibrium for the tracer particle 
when the mass ratio becomes small, Eq. (|1(JH . In Fig. 
EJwe show numerically exact solution of the problem for 
cases (1-3). These solutions, obtained using Eq. (J2J for 
finite values of e, show a good agreement between numer- 
ical results and the asymptotic theory. The Levy equi- 
librium for the tracer particle is not sensitive to precise 
shape of the velocity distribution of the gas particle, and 
hence like the Maxwell distribution is stable. 

We now consider a Fokker-Planck equation which de- 



scribes the evolution of the tracer particle PDF W(V m ,t) 
towards the Levy equilibrium Eq. (|10() . The equation is 
of fractional order and is obtained using a small e expan- 
sion of Eq. (JJJ (details to be published) 



dw(y M ,t) 

dt 



A-a/2 Q\Vm\ c 



dVl 



M 



(11) 

In Eq. the Riesz fractional derivative was used 
and the dissipation term is 7 = 2eR. A generalized Ein- 
stein relation 

x/2 



D 



T 



r(l + a) \M 



7, 



(12) 



yields the relation between the transport coefficients D 
and 7. When a = 2 the Einstein relation is recovered. 
Note that [E 0> IS El investigated related fractional pro- 
cesses based on a stochastic approach (e.g. Langevin Eqs. 
with Levy noise) . In those investigations dissipation and 
fluctuations were treated as though they are independent, 
hence the equilibrium obtained there differs from ours. 

Model 2 The question remaining is wether Levy equi- 
librium a general feature, which might be obtained from 
other collision models. Specifically, it is interesting to see 
if Levy equilibrium is compatible with a non-linear Boltz- 
mann equation approach. Since one may suspect that the 
Levy behavior obtained so far is limited to linear Boltz- 
mann models. For this aim we investigated the one di- 
mensional driven inelastic Maxwell model (DIMM) . This 
model was investigated extensively in recent years in the 
context of inelastic gases assuming finite variance bound- 
ary conditions (see details below) 0, 0] . Our goal is to 
investigate DIMM in the quasi elastic limit showing that 
Levy statistics describes the equilibrium, the Maxwell- 
Gauss distribution is recovered in the proper limit. 

First consider the inelastic Maxwell model in the ab- 
sence of external driving forces W(V, t) — I(V, W), where 
the non-linear collision integral is 



I(V,W) 



1 f°° 

-W(V,t) + - 

P 7-oo 



W(u,t)W 



qu 



P 



,t)di 



(13) 

In Eq. (|13fl p = (r + l)/2 and q = 1 — p where r is 
the restitution coefficient < r < 1. The kinetic scheme 
describes a situation where momentum is conserved dur- 
ing collision events, while energy is conserved only when 
r = 1. If r < 1 the steady state solution of Eq. I|13|) 
is W SS {V) = S(V), reflecting the loss of energy during 
collision events. Note that for elastic collisions r = 1, 
any initial velocity distribution is a steady state solu- 
tion. This is expected (and not informative) since two 
identical ID elastic particles, exchange their velocities in 
collision events. 

Let W(k,t) be the FT of W(V,t). The boundary con- 
ditions we will consider are 
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for small k. Since W(V, t) is a non-negative PDF we 
have < a < 2. Using the Boltzmann equation (|13fl it 
is easy to show that in the elastic limit r = 1, ^ = 

0. For the standard case a = 2 considered in Eq. 
(|14f> simply reflects energy conservation, i.e. (V 2 ) is a 
constant of motion. For a < 2, (|F| a ) describes the width 
of the probability packet, which for elastic collision is a 
conserved quantity. 

As mentioned, when r < 1 the inelastic collisions will 
shrink any initial probability packet to be concentrated 
on V = 0. Similar to previous work |lll| an infinites- 
imal heating term is added to the equation of motion, 
which compensates the energy loss. We will consider the 
boundary conditions described in Eq. I|14|) . while ^l| 
considered the Gaussian case a = 2. To obtain behav- 
ior compatible with Eq. (|14|l we consider the fractional 
DIMM 



dW(V,t) 
dt 



d a W{V,t) 



I(V,W). 



(15) 



d\V\ a 

It is more convenient to consider this fractional equation 
in Fourier space, this yields the non-linear and non-local 
equation 

W {k, t) + {l + D a \k\ a ) W [k, t) = W {pk, t) W {qk, t) . 

(16) 

For our aim this equation gives the definition of the frac- 
tional derivative in Eq. (|15|) . Our aim is to investigate 
the steady state solution of this equation in the quasi 
clastic limit when D a — > and r — > 1. This limit is 
taken in such a way that the boundary condition Eq. 
HI 41 is satisfied. 

Using the condition d ^g}°" > = 0, and Eqs. UUUHJ we 
obtain 



D n 



(\v\ a ) 

r(l + a 



{l-p a -q° 



(17) 



Without loss of generality we may set now (|U| a ) = 1. 
For a — 1 we have D a — 0, while for < a < 1 D a 
obtains negative values. The case a — 1 marks the tran- 
sition between a finite {a > 1) and infinite {a < 1) first 
order moment of velocity J_ |V|W(V, t)dV. For a < 1 
no steady state is obtained, since the dissipation due to 
collisions is not strong enough to compensate the heating. 
Our results in what follows are restricted to 1 < a < 2. 
Steady state solution of Eq. (fTtjll satisfy 

(1 + D a \k\ a ) W ss (*) = Wss {pk) Wss {qk) . (18) 

An iteration method is used to obtain the solution, let 
ijj{k) = In [W ss (k)] , and using Eq. (|T%j) we have 

if, (fc) = %j> {pk) + i> {qk) - In [1 + D a \k\ a ] . (19) 

The boundary condition Eq. I|14fl yields ip{k) ~ 
-\k\ a /T{l + a). The solution of Eq. (JUJ is obtained 
using the iteration rule 

^„+l(fc) = i>n {Pk) + <p„ {qk) - In [1 + D a \k\ a ] , (20) 



where \im n - >oa ip n {k) = ip{k) and the 'initial condition' 
is tpo{k) = — In [1 + D a \k\ a ]. Using these rules and some 
algebra involving series expansions, we find 



~ (-l)" +1 \k\ an Dl 
4^ n{\ -q an -p an )' 



(21) 



where the condition 1 < a < 2 was used. Inserting Eq. 
ifTTjl in Eq. (J2U we obtain 



a(l — r) 



\k\ 



r(i 



8 



[T(l + a)Y 



0{\k\ ia ) 



(22) 

The interesting thing to notice, is that the second term on 
the right hand side of Eq. (1221) vanishes when the elastic 
limit r — ► 1 is considered. Inserting Eq. l|17fl m Eq. 
(|21|l one can show that in the elastic limit lim r ^i ip(k) = 
— \k\ a /T (1 + a). Hence the steady state characteristic 
function is a stretched exponential 



lim W ss {k) = exp 

r— »1 



\k\ a 



r(i + o) 



(23) 



the inverse FT of this equation yields the symmetric sta- 
ble Levy density [2| . It is rewarding to find that Maxwell 
and Levy equilibrium are obtained only in the elastic 
limit, thus conservation of energy in the collision events 
is related to the Levy-Maxwell behavior. Far from this 
limit results not directly related to the Gauss-Levy cen- 
tral limit theorem are obtained, Eq. (|21|l . 

To conclude, we demonstrated the relation between 
equilibrium properties of very different types of collision 
models and Levy statistics. Thus stable behavior tran- 
scends details of individual models and hence I suspect 
can be found in other types of collision models. To sup- 
port the idea that Levy velocity distribution might be 
found in other models, we note the interesting work of 
Min et al 0| who used numerical simulations of a long- 
range interacting vortex model, and showed that distri- 
bution of velocity fields are Levy stable. 
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